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Abstract 

I calculate heavy-light decay constants in a nonrelativistic potential model. The 
resulting estimate of heavy quark symmetry breaking conflicts with similar estimates 
from lattice QCD. I show that a semirelativistic potential model eliminates the conflict. 
Using the results of heavy quark effective theory allows me to identify and compensate 
for shortcomings in the model calculations in addition to isolating the source of the 
differences in the two models. The results lead to a rule as to where the nonrelativistic 
quark model gives misleading predictions. 
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1 Introduction 



The nonrelativistic quark model is one of the oldest and most successful models of hadronic physics. This 
success is somewhat puzzling in that it persists even when the model is applied to light quark hadrons, 
where the dynamics are dominantly relativistic. Perhaps more puzzling is that relativistic corrections to 
the nonrelativistic quark model do not to substantially improve the model's predictions for spectra 0. 
Some (but not all) of the ideas of the nonrelativistic quark model for heavy-light systems gain a stronger 
theoretical basis through heavy quark effective theory (HQET). In this work I show how the nonrelativistic 
quark model can be used in conjunction with HQET to calculate heavy-light decay constants. By doing 
the same calculation with a semirelativistic potential model, I show how relativistic extensions of the simple 
quark model can make a dramatic improvement in some types of calculations. This, in turn, indicates which 
nonrelativistic quark model calculations should not be trusted. 

Before turning to the model calculations it is important to understand what HQET tells us about decay 
constants, since HQET provides the only results that follow directly from QCD. The application of the 
ideas of HQET to heavy-light decay constants preceded the development of the effective theory itself. The 
nonrelativistic quark model led to the prediction that heavy-light decay constants follow the scaling behavior 
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Later, Shifman and Voloshin || and, separately, Politzer and Wise calculated the leading-logarithmic 
corrections to Eq. ([!]) 
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in a model-independent manner, i.e., following directly from QCD in the limit where the heavy quark mass 
goes to infinity while the the QCD scale remains fixed. 

The above relation is of both theoretical and practical interest. Theoretically, Eq. (||) is interesting 
because it is a model-independent prediction of QCD in a well-defined limit. Practically, it is interesting 
because fa is an input to other calculations, such as B°-B° mixing. Unfortunately, a direct measurement 
of fs through leptonic decay will be extremely challenging because of the very small branching ratio and 
difficult signature. A measurement of fo, on the other hand, is much more feasible. In fact, measurements 
°f fa,-, which is related to fa by flavor SU(3), are already available 0,0,0, albeit with large errors. 

Unfortunately, in the real world the bottom and, particularly, the charm quark masses are quite finite 
compared to the QCD scale. It is therefore necessary to consider the finite-mass corrections to Eq. (|2|). The 
predictive power of the effective theory vanishes when the leading-order finite mass corrections the decay 
constants are included. This means the size of the corrections must be estimated using lattice QCD or some 
model. This problem has been studied extensively on the lattice where results indicate that the 

corrections to the heavy quark limit for fa are O (€/%), which corresponds to a subleading heavy quark term 
of size (1 GeV)/mQ. QCD sum rules are consistent with these lattice results. The large correction is 
surprising when compared to what one would naively expect from the nonrelativistic quark model, something 
like (0.3 GeV)/mQ. Naive estimates can miss factors of three, of course. It is necessary to do an explicit 
calculation to see that the nonrelativistic quark model really conflicts with the lattice calculations. 

This work uses two simple potential models to explicitly calculate decay constants in the heavy quark 
limit and beyond. The first, hereafter referred to as the "nonrelativistic quark model" is based on the 
Hamiltonian of the Isgur-Scora-Grinstein-Wise (ISGW) model [[l3) in the heavy quark limit. This model is 
very simple, in contrast to lattice methods, which are rigorous, but also exceedingly complicated. Even if 
the lattice is able to provide precise answers to the structure of hadrons, it is useful to find simple pictures 
which describe the important physics. The second model, the "semirelativistic quark model" is an simple 
generalization of the first. The difference is substitution of the relativistic form for the kinetic energy for the 
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nonrelativistic form used in the ISGW model. This simple change involves subtleties which are discussed in 
the body of the text. 

One might reasonably ask, given a willingness to use these models for calculations, why bother with 
HQET at all? There are several reasons. The first is that HQET provides some checks on the calculation. 
At subleading order, HQET does show that there is a term missing from the model calculation. Fortunately, 
it is a term that may be added by hand. A second reason is that HQET allows the inclusion of radiative 
corrections in a rigorous manner. Finally, and perhaps most importantly, HQET provides a detailed way 
to compare models. When models differ in their predictions, it is desirable to isolate the regions in which 
they differ. Unfortunately, when one takes apart two different models to compare, it becomes a matter of 
comparing apples and oranges. By calculating nonperturbative matrix elements that arise in HQET, the 
two models can be compared in a physically meaningful way. 

The next section reviews the HQET predictions for decay constants to subleading order in I/tuq. The 
following sections describe the calculations in the nonrelativistic and semirelativistic models. I then compare 
the results of the two models and discuss the implications for other nonrelativistic quark model calculations. 
The appendix describes the numerical methods I used to do the calculations. 



2 HQET for meson decay constants 

The heavy quark effective Lagrangian |14|,|f5|], 



AtQ = KD ■ vh v , 



(3) 



is by now well known. For a review which includes an extensive discussion of decay constants, see Ref. flT6f . 
The spin and heavy quark mass symmetries of the heavy quark limit are manifested by the lack of gamma 
matrices and masses in Eq. (|J). The usual definition of the pseudoscalar decay constant, Jm, of a Qq meson 
M with four-momentum p is 

(0|A M |M(p)) = if M p„, (4) 

where is the axial current. Throughout this work M (M*) represents a heavy-light pseudoscalar (vector) 
meson with a heavy quark Q and a light antiquark q. Using the symmetries of Eq. (|^), one can see that in 
the heavy quark limit 

/mV^m = F, (5) 

where F is a universal dimensionful parameter of QCD. This parameter depends on the nonperturbative 
sector of QCD, so it is not currently calculable from first principles, (ft is calculable on the lattice in 
principle.) In the symmetry limit, i.e., when the bottom and charm quarks are taken to be infinitely 
massive, the decay constants of the D, D* , B and B* are determined by F. 

The discussion so far ignores radiative corrections. When the leading-logarithmic radiative corrections 
to the axial current in the heavy quark limit are included, the result becomes 



fM\fm M 



(6) 



where (3q — (33 — 2n/)/3. The general form of the result is C(n)F(fi), where C(/i) is the perturbative 
coefficient to the low-energy parameter Since physics does not depend on the choice of scale, the 

/i-dependence of the product must vanish. 

At subleading order in 1/toq and including leading-log radiative corrections, the Lagrangian grows [ p~7| , 
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A third term at order 1/toq whose matrix elements vanish due to the equations of motion has been omitted. 
The first correction term is the leading part of the kinetic energy of the heavy quark. Its perturbative 
coefficient is unity because of reparameterization invariance p(i(| . The second correction term arises from the 
heavy quark's non-zero chromomagnetic moment. These terms give rise to corrections to the decay constant 
through modifications of the meson wave function and of the heavy-light current. When these effects are 
included, the simple result in Eq. (||) becomes (ignoring radiative corrections for simplicity) |21j| 



/ma/ to m = F < 1 
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where dj\t = +3 (—1) for pseudoscalar (vector) mesons. Here the effect of the modification to the meson 
wave function due to the kinetic energy term and the chromomagnetic term are parameterized by G\ and 
G2, respectively. The finite difference between the heavy quark momentum and the heavy meson momentum 
gives rise to the final term. G\, Gi and A are dimensionful parameters of QCD which, like F, cannot be 
calculated in perturbation theory. 

Unlike G\ and Gi, the parameter A is directly related to other heavy quark processes. The difference 
between the mass Tom of a heavy-light meson and the corresponding heavy quark mass tuq is conventionally 
defined as 

A = m M - rriQ. (9) 

A is related to A by |l§ 

A = A-m 3 , (10) 

where m q is the light quark mass. Since the current masses of the up and down quarks (5-10 MeV) are 
considerably smaller than estimates of A (typically 300-700 MeV), A is usually taken to be equal to A. A 
subtlety which arises in the model calculations in the following sections makes this distinction important. 
Including the leading- log radiative corrections to Eq. (||), 
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An important point about the above complicated expression is that the Gi(/i) term has the same pertur- 
bative coefficient as F(fi) because of reparameterization invariance, while the Gi((a) term gets a non-trivial 
perturbative coefficient. 

At this level the results of HQET have lost their beautiful simplicity. Unfortunately, they have also lost 
their predictive power because the four decay constants f D M, f B (*) are given in terms of four unknown 
parameters F, Gi, Gi and A. Of these, only A is obtainable from other heavy- meson processes in principle. 
The task for the model calculations is to estimate these parameters. 



3 Nonrelativistic Model Calculation 



In the nonrelativistic quark model, meson decay constants are given by (22,23, 24|,p5|. 



/a/ v / ^= V12|^M(r = 0)|. 

The "non-relativistic quark model" for this paper is a heavy constituent quark 
antiquark q obeying the Hamiltonian 
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Figure 1: Nonrelativistic wave function. 



The last three terms in the Hamiltonian are the quark-antiquark potential. The first and third represent the 
coulomb-like and hyperfine effects of single gluon exchange, respectively. The linear term is a phenomeno- 
logical spin-independent confining potential. A more general Hamiltonian would also include spin-orbit 
coupling terms. I have omitted such terms because all of the calculations in this work involve only S- wave 
states for which spin-orbit contributions vanish. Without the hyperfine term, the Hamiltonian is that of the 
ISGW model fll3||. W ith the hyperfine term, the model is closely related to the updated model of fsgur and 
Scora (ISGW2)[§7). It should be noted that this model differs from ISGW and ISGW2 in that I use exact 
(numerical) wave functions, while the others use simple variational wave functions. While variational wave 
functions are useful for the calculations in ISGW and ISGW2, which involve overlaps of wave functions, they 
are not appropriate for decay constants, which are sensitive to the wave function at a single point. 

As the heavy quark mass is taken to infinity, expectation values of p and S 3 (r) remain of order of the 
QCD scale. In this limit the above Hamiltonian reduces to p8[ 
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Solving the Schrodinger equation 



for the ground state wave function and comparing with Eq. ( |12| ) gives 

F=x/l2>NR( r ^ 0) |. 



(14) 
(15) 
(16) 

This model calculation explicitly obeys the mass and spin symmetries of the heavy quark effective theory. 
Solving Eq. @ numerically gives F — 0.55 GeV 3 / 2 . I have used the parameter values m q = 330 MeV, 
a s — 0.5, a — 0.18 GeV 2 from Ref. p3| . Figure [l] displays the calculated wave function. Unfortunately, 
numerical calculations such as this tend to obscure the dependence of the results on the input parameters. 
This is particularly important when trying to establish agreement or disagreement between different types 
of calculations. Figure |^ makes the parameter dependence of the result more explicit by displaying the 
dependence of F on the input parameters within ±50% of each nominal value. 
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Figure 2: Parameter dependence of F calculated using the nonrelativistic model. The nominal values (.To) 
are (solid line) a s — 0.5, (dashed line) m q = 0.33 GeV, and (dotted line) a = 0.18 GeV 2 . 



In order to calculate the decay constant to subleading order in 1/toq, it is necessary to reintroduce the 
heavy quark kinetic energy and hyperfine interactions to the Hamiltonian. The wave function can be written 
in an expansion in powers of 1 /rag as follows 
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The functions and </>jjp arise from the effects of the kinetic energy and spin-spin hyperfine terms, 

respectively. They have been defined to be independent of rriQ. Note that A(S C • Sd) = (—3, 1) for (M,M*) 
is —cLm, which was defined in the previous section. Simply solving the Schrodinger equation including the 
1/toq terms leads to contributions of the subleading-mass terms to all orders in I/tuq. The functions 
and 0jjf can ^ e isolated using perturbation theory, where 
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In Eqs. ( |18| ) and (19) the set of functions (f>oo, 4>i, 4>2- ■ ■ represents the complete set of solutions to Eq. (15) 
with the unperturbed Hamiltonian in Eq. (|l4|). In practice it proves easier to numerically find the piece of 
the solution to the full Hamiltonian linear in I/tuq than to directly apply Eqs. (18) and (|l9|). 
The constants G\ and Gi defined in Eq. (0) are related to the wave function corrections by 
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The quark model calculation to order 1/mg reproduces the form of the heavy quark result in Eq. (||) with 
the exception of the term proportional to A, which is absent in the model calculation. This missing term 
manifests one of the limitations of the constituent quark model. It can be understood as follows: the factor 
A = A — m q arises in a process with q 2 = m 2 M , which is large compared to Aq^fj . The relevant light quark 
mass m q should therefore be the current quark mass, which means 

A = A - m c q n " cnt ~ A. (22) 

The quark model only knows about constituent quarks, however, which would give 

quark model _ ^constituent g (23) 

This facet of the quark model calculation is wrong. Fortunately, the deficiency can be compensated for by 
manually including the A term. 

A more serious problem arises in the calculation of the hyperfine correction to the wave function (<^Sp)- 
A straightforward evaluation of the sum in Eq. ( |l9| ) shows that </>jjp (and consequently G2) diverges. The 
delta-function potential is too singular for the Schrodinger equation so the wave function at the origin 
diverges, even at leading order in perturbation theory. Although it is possible to regulate this singularity 
through a variety of methods, the resulting calculation depends critically on the method chosen. Since the 
effect of the perturbation on the wave function (eigenfunction) is infinite, one might naively expect that the 
effect of the perturbation on the mass (eigenvalue) would also be infinite. Then the correction to the heavy 
particle mass, which is measurable through the B-B* mass splitting, could be calculated in the regularized 
theory and subsequently be used to fix the regularization parameter by fitting to the measured mass splitting. 
Unfortunately, this procedure fails because the effect of the hyperfine perturbation gives a finite eigenvalue 
correction, even though it gives an infinite eigenfunction correction. 

I will assume that the hyperfine contribution, and subsequently G2, is negligible compared to the kinetic 



energy contribution with the following justifications: Qualitatively, one can compare the terms in the Eq. (19) 
sum with the terms in the Eq. Jl§| ) sum. In both cases the first terms in the series, i.e., the contributions of 
the lowest-lying excited states, are larger in magnitude than any other terms. Comparing only these first few 
terms, the kinetic energy perturbation is much larger than the hyperfine perturbation. However, the large-n 
terms in the hyperfine sum fall only like 1/n, so the sum diverges, whereas the terms in the kinetic sum 
fall quickly enough for the sum to converge. From this it seems plausible that an appropriately regularized 
calculation will yield | G2 1 < \Gi\- Furthermore, two QCD sum rule calculations [^l],[l^] give | G2 1 "C \G\\. 
Therefore, it is reasonable to assume that the hyperfine interaction can be neglected in this calculation. 

Fortunately, the dominant G\ term is easily calculable in the nonrelativistic model. The effect of (re- 
introducing the heavy quark kinetic energy can be incorporated in the usual way by substituting the reduced 
mass m re( j, 

1 1 I 

(24) 



m red m q m Q 

for the light quark mass m q . However, simply making the substitution introduces corrections to all orders in 
1/uiq. This is a problem because the heavy quark kinetic energy p 2 /(2itiq) is only correct to leading order 
in l/mQ.|^| Taylor-expanding 4>(0) as a function of m rec j 
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T Which, of course, leads one to wonder about higher-order corrections to the light quark kinetic energy, 
which do not converge, since typical values of p are of the same order as the constituent light quark mass 
m q . Concerns such as these inevitably lead to a model with relativistic light-quark kinematics such as the 
one in the following section. 
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Figure 3: Parameter dependence of G\ calculated using the nonrelativistic model. The parameters (x) are 
a s (solid line), m q (dashed line) and a (dotted line). 



yields the following expression for G\ 



m l dcfroo (r) 



(26) 



0oo(O) dm q 

Numerically, it is easier to treat the heavy quark kinetic energy as a perturbation, as described above. 

The numerical calculation gives G\ = — 0.14 GeV. Fig. || displays the parameter dependence of the 
calculation, showing that varying the parameters does not allow for values of G\ much larger (in absolute 
value) than 0.3 GeV. 



4 Semirelativistic Model Calculation 

The idea of the semirelativistic model is to remove the most obviously incorrect part of the nonrelativistic 
quark model, the nonrelativistic form of the kinetic energy for the light quark. Take the nonrelativistic quark 
model Hamiltonian in the heavy-quark limit, Eq. (|l4|), and make the substitution 

2 

V + (27) 



2m q 



The resulting wave equation, 

p 2 + m 2 ^ = „ ^ ( 28 ) 

is known as the spinlcss Salpctcr equation []29j]. It follows from the full Bethe-Salpeter equation in the 
spin-independent and instantaneous-interaction approximation. The spin-independence is justified by the 
heavy-quark limit. The instantaneous-interaction approximation is a limitation of the model. Duncan, 
Eichten and Thacker have shown ]30] | that the spinless Salpeter equation produces wave functions which are 
very similar to those obtained from lattice calculations. 

If the substitution of the relativistic kinetic energy is the only change made to the model of the previous 
section, the potential in Eq. (j|^) is 

V = + ar. (29) 

3r 
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Unfortunately, the resulting solution to Eq. (|28j) diverges at the spatial origin, which results in an infinite 
value for F when calculated with Eq. ([l6]). One might be tempted to ascribe this divergence to the phe- 
nomcnological part of the potential. However, the divergence depends only on the coulombic part of the 
potential; it is independent of the phenomenological linear term. Wave function divergence at the spatial 
origin is actually a general problem affecting relativistic wave equations. For example, the solution to the 
Dirac equation for the Coulomb potential behaves like 

~ {2mar) Vr ^- 1 (30) 



for small r. While the divergence in Eq. (|30|) is very weak, the divergence of the solution to the spinless 
Salpeter equation with the potential in Eq. (^9|) is much stronger, 

*^ r -^r ; (31) 

as can be seen with the methods of Ref. j3]J . 

The singularity in the wave function is clearly related to the singularity of the 1/r potential. If we 
consider instead the one- loop single gluon exchange potential |3^,[33|], the net effect is to replace the constant 
value of a s in Eq. (|l|) with the one- loop running value of a s (X/Ar). Leaving the phenomenological linear 
term unchanged, the potential is now 

v = _^iMM +art (32) 

where 

4-7T 

as(1/Ar) = fl,ln(l/AV) ' ( 33 ) 

which has a much milder singularity at the origin. The resulting wave function still diverges, but only 
logarithmically. Again following a derivation similar to that in Ref. |3l}| , one can show that the small-r 
behavior of the solution to the spinless Salpeter equation with the potential in Eq. ( |32|) is 

^(r^O)~[-ln(Ar)] 4/3flo . (34) 

The physical decay constant is a product of a perturbative coefficient which depends on a scale fi with 
the low-energy parameter F(p), as is shown in Eq. (^[), which is correct to leading-log order. The logarithmic 
behavior of the wave function in Eq. ( |34| ) is of the right form to cancel the ln(/i) dependence of the perturbative 
coefficient that would be obtained if we had only considered single gluon exchange (i.e., the vertex correction) 
in the perturbative coefficient in Eq. @. This is as it should be, since the solution to the wave equation 
can be considered an infinite series of single-gluon exchanges. The full one-loop perturbative calculation also 
includes the propagator corrections for the light and heavy quarks, but these effects are not present in this 
model. A better model would produce the full ln(yu) dependence of F(fi). 

In the present case, the correct quantity to compare with the nonrelativistic wave function at the origin 
is the semirelativistic wave function at the origin without the logarithmic divergence, which should cancel 
with the fj, dependence of the perturbative correction in the full calculation. Then the quantity 

4> SR (r) 

F ° = ^oln(l/Ar)4/3/3o ( 35 ) 

should be compared to F as calculated in the nonrelativistic model. 

A subtlety arises in the one-loop potential because a s (l/Ar) diverges for r ~ A -1 . This unphysical 
behavior arises from the nonperturbative nature of QCD at long distances, and as such should be swept 
into the phenomenological part of the potential. I have followed the procedure used by Peskin and Strassler 
|34|| to smoothly turn off the running of a s at long distances. The prescription is to make the substitution 
Ar — * Ktax\h(Ar / k) in the running of ct s . The results are insensitive to the precise value of the parameter k, 
which I set to 0.5 for the results I present here. 
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Figure 4: Semirelativistic wave function calculated using 10 (dotted line), 15 (dashed line) and 20 (solid 
line) pseudohydrogenic basis functions. The inset shows that the numerical calculation fails to converge at 
the origin, where the wave function diverges logarithmically. 



Fig. ^ shows the numerical solution to the spinless Salpeter equation with the one-loop potential. The 
values of m q and a are the same as in the previous section. There is a subtlety in choosing m q in this 
model. In one picture the constituent quark mass arises from the relativistic "jiggle" of the light quark in 
the hadron. In another picture, the constituent quark mass arises from chiral symmetry breaking. Although 
these schemes are not necessarily mutually exclusive, the former requires using the current light quark mass 
in this model, while the latter requires using the constituent mass. Here I have chosen the latter option. It 
should be noted, however, that the results do not depend very strongly on the light quark mass, so choosing 
the former option would not qualitatively change the results. I have chosen A = 237 MeV so that the 
one-loop potential (Eq. (|3^)) is the same as the orig inal potential (Eq. @) at r = 1 GeV -1 . Fig. | displays 
the sensitivity of resulting value of Fq to the model parameters. The central value, Fq = 0.67 GeV, is only 
about 20% higher than the value of F obtained in the nonrelativistic model. 

Calculating the subleading terms A, G\ and G2 is quite similar to the nonrelativistic calculation. The 
A term has to be included by hand, just as before. The hyperfine effect, G2, which was problematic in the 
nonrelativistic model, is also problematic in the semirelativistic model. Even though the one-loop potential, 
with its milder r — > singularity, helped deal with the divergence at the origin of it does not alleviate 
the additional singularities in the hyperfine potential. The one-loop hyperfine potential |B2,B3], 
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contains terms just as singular as the tree-level hyperfine potential. This means that even after the original 
divergence at the origin is regulated, the hyperfine potential will introduce a new divergence. This result 
can be anticipated from the perturbative calculation in the effective theory, where G2 (jJ<) gets a perturbative 
coefficient in Eq. ( pT| ) beyond that of F(p). As a rule, terms which have renormalization coefficients with non- 
trivial ln(/i)-dependence diverge in the semirelativistic model. G\, which is protected from renormalization 
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Figure 5: Parameter dependence of Fq calculated using the semirelativistic model. The nominal values (xq) 
are (solid line) A = 0.237, (dashed line) m q = 0.33 GeV, and (dotted line) a = 0.18 GeV 2 . 

by reparameterization invariance, does not diverge in the model calculation. 

Although the semirelativistic model Gi calculation seems to be more tractable than the similar problem 
in the nonrelativistic model, the calculation is extremely sensitive to the small-r dependence of the wave 
function. This is precisely where the numerical method breaks down, so the calculation is not technically 
feasible. Neither the nonrelativistic nor semirelativistic models in present form give definite predictions for 
G2. Fortunately, as argued in the previous section, indications are that Gi is negligibly small compared to 



The kinetic energy term, G±, can easily be calculated by treating the p 2 /2mQ term as a perturbation. 
(The interpretation as a reduced mass effect mentioned for the nonrelativistic model does not translate to 
the semirelativistic model, so Eq. (|2^) no longer holds.) The result is displayed in Fig. 0, once again showing 
dependence on the various parameters. 

Here the two models give dramatically different results. G\ is several times larger in the semirelativistic 
model than it is the nonrelativistic model. Also the parameters are correlated such that the two models 
cannot be made qualitatively similar by changing any combination of the parameters. 



In Table |l] I have summarized the results of the two model calculations, including for pseudoscalar 
mesons, gp. (In general, the pseudoscalar meson's gp will be different from the vector meson's gv due to 
the effects of the hyperfine operator, Gi. Since the model calculations neglect the hyperfine contribution, 
9m = 9p = 9v f or these calculations.) I have used A = m q = 0.33 GeV to calculate gp. Note that gp in 
the nonrelativistic calculation is equal to the naive guess from the introduction. This is actually fortuitous, 
because the contribution of the A term, which accounts for half of the value, is not included in the naive 
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5 Discussion 



It is convenient to define a quantity gM, such that 




(37) 
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Figure 6: Parameter dependence of G\ calculated using the nonrelativistic model. The parameters (x) are 
A (solid line), m q (dashed line) and a (dotted line). 



Table 1: Comparison of the nonrelativistic and semirelativistic models. 



F (F ) [GeV 3 / 2 ] 


Gi [GeV] 


g P [GeV] 


nonrelativistic 0.55 


-0.14 


-0.31 


semirelativistic 0.67 


-0.95 


-1.12 
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model. As stated in the introduction, lattice calculations indicate that gp ~ 1 GeV. The semirelativistic 
calculation is consistent with that result, but the nonrelativistic calculation is not. The difference is due to 
the G\ contribution. 

The nonrelativistic and semirelativistic values of G\ differ not just quantitatively, but qualitatively. This 
qualitative difference would be completely obscured by a model comparison done in the traditional way, i.e., 
by calculating only the decay constant and including the heavy quark mass effects to all orders. The heavy 
quark mass suppresses the effects of the G\ term in the decay constant itself. The (heavy quark suppressed) 
large difference at subleading order also tends to compensate the smaller difference between the two models 
at leading order. 

The origins of the discrepancy between the two models can be understood as follows: For small p, the two 
Hamiltonians are the same. For large p, however, the kinetic energy term grows like p 2 nonrelativistically, but 
only like p relativistically. This means that the semirelativistic Hamiltonian is less confined in momentum 
space than the nonrelativistic Hamiltonian, i.e., the semirelativistic wave function is more spread out in 
momentum space than the nonrelativistic wave function. Because the wave functions are normalized, an 
increase of the wave function at large momentum must be compensated by a decrease at small momentum, 
so the difference between the wave functions tends to cancel for ij){r — 0), which can be written in p-space 
as 

tp(r = 0) = J d 3 p^(p). (38) 
Gi, however, is proportional to 4>KE( r — 0), which can be written as 



W = E it^lL I d V<(p')|^V)- (39) 



The p 2 factor in the integral emphasizes the large-p differences in the wave functions, making G\ a sensitive 
probe of the large-momentum tail of heavy-light wave functions. This, in turn, leads to the following rule: 
Quantities which are sensitive to the large-momentum shape of wave functions are dramatically underesti- 
mated by the nonrelativistic quark model. 

This rule has implications for other processes. In particular, it indicates that nonrelativistic quark model 
calculations of processes at large momentum transfer seriously underestimate the overlap of meson wave 
functions. An important example that has received much interest lately is the process B — > K*-f. In the 
B meson's rest frame the K* has rj 1,3 GeV of momentum, which is large compared to the typical widths 
of meson wave functions in the nonrelativistic quark model. This means that the overlap is dominated by 
the tails of the wave functions, which I have just shown to be poorly described by the nonrelativistic quark 
model. 

This work not only provides an explanation for the conflict between the nonrelativistic quark model and 
other estimates for the heavy-quark symmetry breaking behavior in decay constants, it also suggests a qual- 
itative solution to an earlier conflict: In Ref. |p8| , I calculated heavy-quark symmetry- violating corrections 
to form factors in B — ► D^*Hu transitions. The predictions for the effects of the heavy-quark kinetic energy 
operator in were an order of magnitude smaller than a QCD sum rule calculation |35) of the same effect. 
This work shows that the nonrelativistic quark model dramatically underestimates the effect of the kinetic 
energy operator, in this case by a factor of 6. In the meantime, Neubert has derived a theorem showing 
the sum rule used in Ref. |3f| overestimate the effects of the same operator. While an explicit calculation is 
needed for both, it appears that the two different types of models should now be in qualitative agreement. 



6 Conclusions 

The nonrelativistic quark model provides a very simple picture of hadronic physics. While the picture 
is clearly too simple, it does yield insight into the structure of hadrons. Calculating decay constants in 
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the nonrelativistic quark model and a simple scmirelativistic generalization shows how the nonrelativistic 
quark model can work reasonably well overall, yet fail to describe important details. This calculation shows 
that the nonrelativistic quark model does conflict with lattice and QCD sum rule predictions for the size 
of heavy quark symmetry-breaking effects in heavy-light decay constants. However, this conflict can be 
removed by going to a similar model with relativistic light quark dynamics. Comparing the two models 
shows that the nonrelativistic quark model should be expected to fail for calculations which are sensitive to 
the large- momentum tails of wave functions. 
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Appendix 



This appendix describes the method I used to obtain the numerical results in the text. While the Schrodinger 
equation can easily be solved with a wide variety of numerical techniques, the spinless Salpeter equation is 
much more difficult . The position-space representation of the spinless Salpeter equation contains the prob- 
lematic V 2 + to 2 , operator, while the momentum-space representation contains a complicated convolution 
integral from the potential. 

These problems can be avoided by using the Rayleigh-Ritz-Galerkin (RRG) method |3^| , which easily 
handles both the Schrodinger and the spinless Salpeter equations. RRG is an extension of the elementary 
variational method. In the variational method, one chooses a state parameterized by A, then minimizes 

£var = (A|W|A) (40) 

with respect to A. For reasonable choices of |A), E^ 1 and |A mm ) form good approximations to the eigenvalue 
and eigenket, respectively. In the RRG method, one chooses an orthogonal set of n vectors, |A, i) i = 1 . . .n, 
then minimizes 

Errg = (*|«|*>, (41) 

where = c,;|A, i). One can calculate wave functions and energy eigenvalues arbitrarily well by choosing 
sufficiently large n. The problem is reduced to the numerically straightforward problems of calculating 
integrals and solving a matrix equation for the Cj's. 

The difficulties with the spinless Salpeter equation can be avoided in the RRG method by breaking up 
the expectation value of the Hamiltonian into kinetic and potential pieces, 

H = T + V, (42) 

then writing Eq. (|4l|) as follows 

^RRG - (*|p)(p|T|p)(p|*) + (*|a:)(a;|V|a;)<a;|*). (43) 

The resulting integrals are straightforward as long as the representations of the |A, i)'s are known in both 
position and momentum space. 

For the calculations in this work I used two different bases: the harmonic oscillator basis and the confined 
pseudohydrogenic basis. The former are standard; the latter were developed in Ref. Q and first used for 
the spinless Salpeter equation in Ref. p9|| . All convergent results are independent of basis. The two different 
bases act as a cross-check. Since the spinless Salpeter wave function diverges at the origin, the two methods 
do not agree in a small region around the origin. The inset in Fig. ^ shows the failure to converge in the 
pseudohydrogenic basis. The same plot with the harmonic oscillator basis is different in the vicinity of the 
origin. Nevertheless, the limiting procedure in Eq. (]3^) provides a finite quantity which is basis- independent. 
All the results in the text are independent of basis. 
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